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Exotic states in the charmonium family are systematically treated in the framework of simplest 
model with an effective coulomb-like interaction of heavy quark and antiquark in the presence of 
static excitation of quark-gluon modes responsible for a nonperturbative term of potential, which 
provides with the confinement of quarks, in terms of bag over the threshold in the excitation spectrum 
of vacuum fields. Once the spectrum has got quite a wide mass gap, it allows us to approximate 
the bag contribution into the potential by a constant value of bag mass at low distances less than 
the bag size. The bag mass can be evaluated in a constituent model. The analysis is given for 
the bag contribution into the distribution over the invariant mass of two pions in the hadronic 
transition between the S'-wave states of bag-quarkonium and heavy quarkonium, that leads to the 
anomaly violating the chiral limit in the region of low invariant masses, which agrees with the 
observational data. Leptonic constants of vector states are investigated in the presence of exotic 
states in the framework of quasilocal sum rules. The extra states allow us to improve the consistency 
of describing the measured widths of leptonic decays for the complete set of vector states in the 
charmonium family. 



PACS numbers: 14.40.Gx, 12.39.Jh, 12.39.Pn, 13.20.Gd 



I. INTRODUCTION 



In the very beginning of charmonium era, it was quite 
spectacular that the potential framework was not only 
qualitatively applicable to the system of heavy quark and 
antiquark, but also quantitatively successful in describ- 
ing the mass spectrum of know members of family as well 
as relevant static characteristics determining mechanisms 
of producing the bound states and their decays. The 
charmed quark is assigned to be heavy since its mass is 
much greater than the scale characterizing the strong in- 
teraction in the framework of Quantum Chromodynam- 
ics (QCD), niQ 3> Aqcd- The nonrelativistic motion of 
heavy quarks with a velocity vq <C 1 provides with the 
multipole expansion in QCD [l], 0, [H , so that in the lead- 
ing order the static potential V(r) can be introduced. At 
short distances between the quarks, the potential can be 
calculated in the perturbation theory 0, H, H, 0] , while 
at large distances it can be approximated by the term 
giving the quark confinement, i.e. the quark bounding 
inside colorless states. We can mention the simplest ex- 
ample of such the potential, the Cornell model || . In this 
model the quarks in color-singlet state are attracted due 
to the coulomb-like interaction with an effective constant 
of "single-gluon exchange" a e s s and confined by the term 
linearly raising with the distance increase. The slope of 
energy increase a is consistent with the linear Regge tra- 



jectories in the spectroscopy of hadrons composed of light 
quarks, 
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where the factor of | is caused by projecting to the color- 
singlet state of quark-antiquark system. The value of 
effective constant has been determined empirically by the 
excitation spectrum in the charmonium family, while a w 
0.18 GeV 2 . Average sizes (r) of charmonium states are 
positioned in the intermediate region of transition from 
the perturbative regime to the nonperturbative one, so 
that the potential is quite accurately approximated by 
the logarithmic function versus the distance [9( with a 
dimensional parameter T determining an average kinetic 
energy in the heavy quarkonium in accordance with the 
virial theorem 1 : 



V((r)) 



(r) 



(2) 



whereas T w 0.38 GeV, so that T ~ Aqcd- 

Therefore, a characteristic square of heavy quark mo- 
mentum in such the quarkonium is determined by the 
product of kinetic energy by the reduced mass equal to 
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A weak dependence of average kinetic energy on the heavy quark 
mass is phenomenologically taken into account in the Martin 
model of potential |10(I . functionally behaving as the power law r k 
at k — » 0, that corresponds to a small perturbation of logarithmic 
form in the region of approximation. 
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mQ/2, i.e. p 2 Ri TOqT, hence, the square of heavy quark 
velocity is equal to 



T 



(3) 



so it is suppressed in the case of heavy quarks. For exam- 
ple, the charmed quark with the mass of m c ~ 1.6 GeV 
has got 

v 2 ~ 0.25, 

and the potential model seems to be quite realistic, since 
relativistic corrections (in the color-singlet state) have 
got a magnitude of Vq. 

The further progress in the development of potential 
approach was related with the introduction of running 
coupling in QCD at small distances in 1-, 2- and 3-loop 
order: the potentials by Richardson [ll[, Buchmiiller and 
Tye [l2j and Kiselev, Kovalsky and Onishchenko 
correspondingly, consistently with linearly raising term 
confining the quarks. This way has allowed us to remove 
a discrepancy between direct measuring the coupling of 
a s and its value, following from fitting the real spectrum 
of quarkonia in the potential model with the running cou- 
pling constant only at the 3-loop approximation 2 . 

However, the theoretical eligibility of potential descrip- 
tion as the framework seemed to become questionable, 
when one clarified that the QCD vacuum has got a 
nontrivial structure, namely, it is composed of quark- 
gluon condensates with the energy scale of the order 
of Aq CD ~ 0.3 GeV Therefore, fluctuations of 

vacuum fields are characterized by the time interval of 
T qcd ~ 1/Aqcd, while the relaxation time of quark- 
gluon fields in the potential description is determined 
by the ratio of heavy quarkonium size to the velocity 
of heavy quark motion, 
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so that by taking into account 
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we get 
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The vacuum fluctuations are irrelevant to the potential 
framework, if the quarks can be considered as static and 



2 The value of coupling constant is inherently related with the 
splitting between the 15 and 25 levels in the quarkonia. At one- 
and two-loop orders the coupling constant from the potential 
model has got a too high value beyond the interval obtained by 
the direct measurement of a s . 



their interaction is practically instantaneous in compari- 
son with the influence of nonperturbative effects, i.e. 



tqcd > tq & 
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In other words 3 , the potential description can be consid- 
ered as rather reasonable and theoretically sound, if only 
the kinetic energy of heavy nonrelativistic quarks is much 
greater the energy scale fixing the quark-gluon conden- 
sates. This fact was recognized by M.B.Voloshin 4 [JlH, 
who gave the estimate for a minimal acceptable value 
of heavy quark mass, for which the description can be 
done in terms of potential models with account of QCD 
condensates: rn m i n ~ 20 GeV, which is in evident con- 
tradiction with the applicability of potential models for 
the both charmed and bottom quarks, since their kinetic 
energy has the same order of magnitude with the scale 
Aqcd, as we have seen above. Thus, the motion of such 
heavy quarks could lead to the excitation of nonzero vac- 
uum fields from the condensates, so that the interaction 
in the system evidently becomes nonstatic, which con- 
tradicts with the success of potential framework in the 
spectroscopy of charmonium and bottomonium. 

The solution of this paradox comes from the quan- 
tum nature of quark-gluon fields. The speculations pre- 
sented above suggest that the spectrum of excitations for 
the vacuum fields is continuous and its threshold is de- 
termined by the energy scale of Aqcd, while the actual 
spectrum could by essentially different. At first, the ex- 
citation spectrum for the vacuum fields can begin with a 
wide mass gap, that is enough in order to substantiate the 
applicability of potential framework to the heavy quarko- 
nium states below the threshold fixed by the mass gap 5 . 
The threshold could be positioned below the threshold for 
the pair production of hadrons with the open charm (the 
production of D-meson pairs in the system of heavy quark 
and antiquark), or it could be very close the continuous 
threshold. Then, the second aspect becomes important: 
the spectrum of excitations for the vacuum fields in the 
presence of heavy quarkonium could be discrete in vicin- 
ity of mass gap, at least in the sense of quasistationary 



3 One could say that the potential would change due to the vac- 
uum fluctuations very slow in comparison with a period of finite 
motion of heavy quarks, i.e. the potential is static, indeed. 

4 In fact, M.B.Voloshin actually formulated the criterium of ap- 
plicability for the multipole approximation in QCD in the 
presence of vacuum fluctuations in the form of condensates: 
tq/tqcd "C 1) while he believed that the necessary constraint 
for the introduction of static potential is the instantaneous in- 
teraction between the nonrelativistic quarks, that oppositely re- 
quires tq/tqcd 2> 1, by his opinion, and hence, it is evidently 
incompatible with the criterium of multipole expansion, So, he 
suggested that the introduction of static nonperturbative poten- 
tial is principally impossible for real charmed and beauty quarks 
except the exotic case of states in vicinity of continuous thresh- 
old. 

5 The potential approach is completely substantiated in the 
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description during the times comparable with the period 
of heavy quark motion in the bound state. Thus, it would 
be quite reasonable to introduce the notion of "bag" for 
the lowest excited state of vacuum fields in the system un- 
der consideration. The bag evidently has got no valence 
degrees of freedom (isospin, charge etc. are equal to zero) 
and it possesses the vacuum quantum numbers (spatial 
and charge parities are positive) except the energy or 
mass. The existence of bag is caused by the presence of 
valence quark and antiquark. In the absence of valence 
quarks, the introduction of bag has no sense. Therefore, 
the bag is not related with notion of glueballs or hybrids 
as particles. In this respect, it would be incorrect to asso- 
ciate the bag with a separate quasistationary hadron-likc 
object, which interacts with other hadrons, for instance, 
with the quarkonium, since the bag itself is not station- 
ary without the quarkonium. Thus, the system of bag- 
quarkonium differs from the hadro-quarkonium recently 
introduced by M.B.Voloshin [17J, ll8|, since the system of 
hadro-quarkonium is considered as the pair of stationary 
objects, which allow the separated existence and interact 
a la the Van Der Waals forces with each other. 

To my opinion, the influence of essentially nonpertur- 
bative phenomenon called the bag on the system of heavy 
quarkonium is reduced to the followings: 

• below the threshold of bag excitation, the heavy 
quarkonium permits the potential description, 
wherein the potential confining the quarks is phe- 
nomenologically approximated by the realistic term 
linearly raising with the distance increase; 

• above the threshold of exciting the bag, the non- 
perturbative contribution is modified: if the heavy 
quarkonium size is less the size of bag Hmg, then the 
bag presence is given by introducing the bag mass 
itself i^bag into the potential, while, if the quarko- 
nium size becomes to exceed the bag size, then the 
linearly raising term of confining potential is acti- 
vated again. 

This situation is schematically shown in Fig. [TJ wherein 
the nonperturbative potential of heavy quarks in the bag 
is pictured. By the way, we assume that further excita- 
tions of vacuum fields are separated from the bag by a 
mass gap, of course. 

The shift of linearly raising asymptote of nonpertur- 
batice term in the presence of bag is determined by a 
interplay between its size and mass. However, it is more 
essential that, if the size of heavy quark-antiquark system 
is less than the bag size, i.e. when one considers the low- 
est states of bag-quarkonium system, then the interaction 
of heavy quark and antiquark "to leading approximation" 
can be written in the form of following potential: 

V cS (r) = ~^ + E hag , (7) 

where af is the effective constant of single gluon ex- 
change between the heavy quarks at the scale character- 
izing the quark-antiquark system. Here we neglect the 
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FIG. 1: The nonperturbative term in the potential energy due 
to the creation of bag in the system of two heavy quark and 
antiquark in the color-singlet state (the bright curve): the 
distance r is measured in fm, the energy has arbitrary units, 
the static limit at zero distance between the quarks is marked 
by label bag, the dashed lines correspond to linearly raising of 
potential at large distances in the presence of bag and in its 
absence. The bag size is "chosen" near 1 fm. 



linearly raising contribution in comparison with the bag 
mass. Thus, we arrive to three-parametric spectral prob- 
lem in the nonrelativistic quantum mechanics: calculate 
the masses of states in the system of bag-quarkonium 
with the coulomb-like attraction at the given mass of 
heavy quark, effective constant a e s s and bag mass -Ebag- 

In Section [TT] we present our treatment of extraordi- 
nary states in the charmonium family in the framework 
of bag-quarkonium system, that allows us to numerically 
estimate phenomenological parameters from the empiri- 
cal data. The parameters are positioned in a region of 
quite reasonable values. In this way we refer to the joint 
table of extra states of charmonia as given in [19j |: Tab. 
HI We exclude the mesons with the hidden strangeness 
Y" s (2175) and beauty Yb as well as the too broad state of 
Y"(4008) presented in the table, and remain them beyond 
our consideration. 

The modification of hadronic transitions from the bag- 
quarkonium states to the ordinary quarkonium with 
emission of pion pair is investigated in Section IIIII The 
presence of bag involves a new term in the amplitude 
in comparison with transitions between the ordinary S- 
wave states of heavy quarkonium, that changes the dis- 
tribution versus the invariant mass of two pions, i.e. it 
can serve as the reason for the observed anomalous be- 
havior of distribution in transitions of ip(3770) — > J/ipmr 
and T(3S I ) — > T(1S)iytt. This fact is the argument in 
favor of our assignment of ^(3770) and T(3S) as basic 
vector states in the system of bag-charmonium and bag- 
bottomonium, respectively. 

Section IIVI is devoted to the analysis of leptonic con- 
stants of vector states in the charmonium family in the 
framework of quasilocal sum rules [201 ] . The study shows 
that widths of observed leptonic decays of vector states 
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state 


mass M (MeV) 


width T (MeV) 


J PC 


decay modes 


production processes 


Ys(2175) 


2175 ± 8 


58 ±26 


1~~ 


0/o(98O) 


e + e~ (ISR), J/tp decay 


X(3872) 


3871.4 ±0.6 


< 2.3 


1++ 


TT + TV~ J/lp^J/lj) 


B -> KX(3872), pp 


X(3875) 


3875.5 ± 1.5 


3.01^ 




D°D°ir° 


B -> KX(3875) 


Z(3940) 


3929 ± 5 


29 ± 10 


2++ 


DD 


77 


X(3940) 


3942 ± 9 


37 ± 17 


J p+ 


DD* 


e+e" -> J/^X(3940) 


F(3940) 


3943 ± 17 


87 ±34 


J p+ 




B -> A"K(3940) 


V(4008) 


4008±^ 


226tE 


I— 


tt + tt~ J/tp 


e+e-(ISR) 


X(4160) 


4156 ± 29 


139^ 3 


J p+ 


D* D* 


e+e" -> J/V'X(4160) 


F(4260) 


4264 ± 12 


83 ± 22 


\— 


ir + ir~ J /tp 


e+e-(ISR ) 


F(4350) 


4361 ± 13 


74 ± 18 


\— 


TT + ir~ip' 


e+e-(ISR) 


Z(4430) 


4433 ± 5 


45l 3 J 


? 




B -» ^±(4430) 


F(4660) 


4664 ± 12 


48 ± 15 


1— 


n + ir~tp' 


e+e-(ISR) 


Y b 


~ 10, 870 


? 


V~ 




e + e~ 



TABLE I: The joint table of cadidates in members of charmonium family: mesons XYZ, as copied from [l£ 



do not contradict with introducting the extra states of 
bag-charmonium. Moreover, the consistency of predic- 
tions with the experimental data is systematically im- 
proved. 

In Conclusion we summarize the obtained results and 
discuss their possible implications to the study of mech- 
anisms for the production and decays of extra heavy- 
quarkonium states. 



II. ANALYSIS OF EXOTICA IN THE 
SPECTRUM OF CHARMONIUM 



The spectrum of energy in the problem with coulomb 
interaction is well known. For the system with the re- 
duced mass m rc d = ^toq, the spin-dependent term due 
to effective single-gluon exchange takes the form 



J 



V S t> = 



2af -L (L ■ S) 1 + \ af\ {3( Sl • n){s 2 ■ n) 



m 



Oi ■ s 2 )} ^3, 



(8) 
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In ([8|) the first term corresponds to spin-spin "contact" 
interaction of quark and antiquark, the second represents 
the spin-orbital interaction, the third gives the tensor 
forces. Here r = T\ — is the relative distance between 
the quarks, the unit vector is n = r/r, L denotes the 



summed orbital momentum of quarks, and S = si ± S2 
gives the summed spin. The masses of levels n +1 Lj 
with n and J being the principal quantum number and 
total momentum are given by 



J 



^4ag-QQK 6+1 £j] = 2m Q + £ b 
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, eff 



off 
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n* \ 3 
\ {L ' S) ~ 4L(£ + l)-3 



(2L±1)((2L± 1) 2 -1) 
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3(L • S) 2 ± I (L ■ S) - L 2 S 2 



r 



whereas 

(L-S) = \{J{J + 1)-L{L + 1)- S(S + l)}, 
L 2 =L(L + l), S 2 = S{S+1). 



The spatial and charge parities, respectively given by 
P = (-l) i+1 and C = (-1) S+L . Particularly, at L = 0, 
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the masses of vector states J PC = 1 are equal to 



^bag-QQ^l] = 2m Q +E b 



m Q f 4 cff 



12ti 3 I 3 5 



: 4n 2 V3 

4 



(10) 



while at L — 1, the masses of spin-triplet J ++ -states are 
equal to 



24n 3 



, cff 




The review of charmonium family is presented in [17l YLty . 
We assign t/;(3770) as the basic extraordinary vector 
state, which is distinguished by the anomalous distribu- 
tion in the spectrum of invariant masses for two pions 
in the transition 0(3770) — > J/ipTtn. One usually con- 
siders that this state is positioned in vicinity of ID level 
of quark-antiquark system cc, and its nonzero width of 
decay into lepton pair is caused by a mixing with the 
nearest ,5- wave level tp' = ip(2S) — 0(3680). Then, one 
can treat it as the ordinary state, and the anomaly men- 
tioned has a reason, which is, for instance, an occasional 
suppression of dominant term in the creation of pion pair 
(see the analysis in [TtJ )■ However, in Section HTT1 we ar- 
gue in favor of that the anomalous distribution versus the 
invariant mass of two pions is the natural consequence, 
if one assigns "0(3770) as the system of bag-charmonium. 
Therefore, we identify 0(3770) with the state 1 in the 
system of bag-charmonium with the quantum numbers 
1 3 Sl 

The other benchmark of exotic state is AT (3872), which 
has got the definite quantum numbers 1 ++ . This state 
has to be identified with the first excitation at the same 
quantum numbers in the system of bag-charmonium, i.e. 
2 3 Pi. Thus, in order to determine the parameters of 
system we have got two levels with the mass difference 
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3m c 
~16~ 



• a 



cff 



17 

36 



• a 



cff 



(12) 



that numerically gives AM = 102 MeV. Eq. (12} allows 
us to write down the effective constant as the function of 
charmed quark mass 



• a 



cff 




1802 AM | 



27 



(13) 



Note that the obtained behavior of effective constant ver- 
sus the charmed quark mass in the limits from 1.3 to 1.8 



GeV is in a good agreement with the 1-loop dependence 
of running coupling constant in QCD 



2tt 



f3 Q ln[fj,/AQCD} 



(14) 



if one puts that the scale /i is determined by the charac- 
teristic momentum of heavy quark in the heavy quarko- 



mum, i.e. 



i 2 = mn T, 



(15) 



wherein T rj 0.4 GeV, A QC d « 190 MeV, and 



Po 



at the number of quark flavors contributing to the renor- 
malization for the scales in problem, rif — 3. The fact 
of consistency between the empirical dependence of (|13|) 
and renormalization group (RG) one of (fT4]) -(|15 p is il- 
lustrated in Fig. [2J wherein the ratio of RG constant 
to the effective one is shown after the normalization at 
77v ~ 1.65 GeV. 



0.9 



1.55 1.6 1.65 1.7 1.75 1.8 

m c , GeV 

FIG. 2: The ratio of running coupling constant in QCD to 
the effective constant, calculated in accordance with (|12|) and 
(|13p : the 1-loop approximation is given by the solid line, the 
2-loop curve is show by dotted line, the dashed line refers to 
the value equal to 1. 

The figure shows that in the wide interval of pole mass 
for the charmed quark, the effective constant coincides 
with the RG running constant within the accuracy bet- 
ter than 5%. The agreement between these two values 
becomes more spectacular in the 2-loop approximation 
for the running coupling constant in QCD 



a 



[2-1] 



(P) 



2tt 



I3 a lii[^/A Q CD} 



Pl Mh/Kqcd] 



where /3i — 51 — 19nf/3, after the appropriate correction 
of Aqcd- However, the usage of 2-loop approximation 
is evidently beyond the accuracy of consideration, since 
the model potential does not take into account the de- 
pendence of effective constant versus the distance in the 
second order over a„. 
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Nevertheless, we can establish that the running cou- 
pling constant of QCD at the scale fixed by the heavy 
quark mass and its kinetic energy slowly depending on 
the mass, allows us to get the empirical value for splitting 
the masses of lightest vector and pseudovector states in 
the system of bag-charmonium. 

A. The pole mass of heavy quark 

It is spectacular (see Fig. [3J that the sum of pole 
masses for charmed quark and antiquark with the bag 
mass is practically independent of the pole mass of 
charmed quark with the accuracy better than 5 MeV, 

£mstab = 2m c + S bag = const., (16) 
after the implication of (|10[) - (fT5|) . 

2m c + -Bbag, GeV 

3,95 , , , , , , 

3.94 - 

3.93 - 

3.92 - 

3.91 - 

3.9 - 

3.89 - 

3.88 - 

1.5 1.55 1.6 1.65 1.7 1.75 1.8 

m c , GeV 

FIG. 3: The sum of bag mass and double mass of charmed 
quark versus the pole mass as the result of fitting the masses 
of lightest vector and pseudovector states of bag-charmonium 
in accordance with (|10[) (|13l) . 

Numerically, we find 

£mstab « 3915 MeV. (17) 

The constant value of this quantity points to that the 
infrared instability of heavy-quark pole mass is exactly 
cancelled by the infrared instability in the static poten- 
tial approximated by the dominant contribution of bag 
mass in the case under consideration. The instability of 
pole mass appears in the perturbation theory of QCD in 
the form of renormalon [111, H3, [H, [H, [H, [2I, [23 , that 
exhibits the factorial growth of coefficients in the pertur- 
bative series as dictated by the presence of infrared pole 
in the coupling constant of QCD at the scale of Aqcd- 
Thus, the observed quantity being the mass of bound 
state for the two heavy quarks and bag is actually sta- 
ble in the infrared region, that we mark by the subscript 
in l]16p . The value of charmed quark mass itself is com- 
pletely determined by the separation of bag mass in the 
approximation considered. 

In this respect it would be interesting to make the fol- 
lowing analogy: at first, consider the light quarkonium 



with the ordinary u- and d-quarks instead of the heavy 
quarkonium. Let us neglect the light quark masses as 
well as the energy of their kinetic motion, if there is no 
orbital rotation. In the basic 5-wave state the mass of 
such the quarkonium with no account of spin-dependent 
interactions is determined by the bag mass 

m(lS) w £ b ag, (18) 

i.e. the state itself is the bag with the valence quan- 
tum numbers given by the light quarks. In the con- 
stituent framework the bag energy fixes the constituent 
mass of light quark by m CO nst. ~ 5-Et>ag- Such the 
bag-quarkonium after the introduction of spin-dependent 
forces is observed as the pseudoscalar ir- and vector p- 
mesons. Therefore, the standard procedure of cancelling 
the contribution of spin-dependent energy leads to 

m(lS) w -{Snip + m w ) » 612 MeV, 

that gives the approximate value of bag mass in accor- 
dance with (fl~8|) . 

Second, the same way can be applied to the considera- 
tion of "constituent model" of quarkonium composed of 
light and heavy quarks. Then, in the basic S'-wave state 
with no account for the spin-dependent forces, such the 
heavy-light bag-quarkonium has got the mass equal to 

M q q[lS] w 771q + TOconst. ~ m Q + - E hag . (19) 

Particularly, for the basic level of D-mesons we get 

Mn[lS] « i £i Rstab « 1958 MeV, (20) 

that should be compared with the experimental value av- 
eraged over the spins for neutral D-mesons, for instance, 

M^ P [1S] = - (3m D . + m D ) » 1972 MeV. (21) 

Comparing (f2Tj) with (|20|) shows that the speculations in 
the framework of constituent model or bag-quarkonium 
allow us to make the calculations with the accuracy about 
15 MeV, characteristic for the potential models in gen- 
eral. 

Thus, in the framework under study we prefer for pre- 
scribing the pole mass of charmed quark to the value 
equal to 

m c Pa \ [£iR St ab - m{lS)] « 1650 ± 15 MeV, (22) 

and analogously, fixing the pole mass of bottom quark 
by the value equal to 

nib ~ m B xp [IS 1 ] — i 777,(15) 

« m c + m c * p [lS] - m c * p [lS] (23) 
w 4995 ± 15 MeV, 
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which surprisingly agree with the values obtained in the 
framework of analysis performed for charmonium and 
bottomonium described by the dominant coulomb inter- 
action with account for corrections both in the pertur- 
bation theor y a nd due to the quark-gluon condensates 
[H, [H, 0, I1H. Then, after the rescaling the effec- 
tive constant of coulomb-like interaction by the renor- 
malization group law to the system of bottomonium, 
we can calculate the mass of basic vector state of bag- 
bottomonium 6 , as shown in Fig. 0J 



M b b °_[l d 5i] w 10355 ± 30 MeV. 



(24) 



Note, that the same result is achieved, if we vary the 
pole mass of charmed quark in the limits from 1.5 to 1.7 
GeV with the consistent change of 6-quark pole mass as 
well as the value of effective constant for the coulomb- like 
interaction. 

M£* [l 3 5i], GeV 

10.5 




5.04 

rrib, GeV 

FIG. 4: The calculation of mass for the basic vector state 
in the system of bag-bottomonium. The vertical band shows 
the permitted values of fa-quark pole mass, the horizontal line 
marks the experimental value of T(3S) mass, the sloped band 
restricts the prediction with the variation of effective constant 
in limits 0.423 < | aj* < 0.484 in the bag-bottomonium sys- 
tem, the bright sloped line gives the prediction for the mass 
at the constant obtained by the rescaling of its value from the 
value for the bag-charmonium in accordance with the renor- 
malization group law. 

At the pole masses fixed above, the coupling constants 
are equal to 



a° ff |cc 



[cc] « 0.479, af[bb]« 0.347 



(25) 



Once such the basic vector state coincides with the position of 
T(3S), possessing the anomalous properties in two-pion transi- 
tion into the low-lying vector levels of bottomonium, while its 
first excitation is positioned at T(4S) with the prediction accu- 
racy of 30 MeV, though the splitting between the further exci- 
tations becomes less than the uncertainty of prediction, so that 
taking into account the width of T(45) comparable with the un- 
certainty of prediction, we do not make certain statement con- 
cerning for the definite correspondence of any exotic state to the 
accepted notation of T(4S). 



Characteristic sizes of heavy quark and antiquark system 
are determined by the "Bohr radius" 



(26) 



2mqtf 



off ' 



which is, for example, approximately equal to a w 0.38 
fm at m c = 1.65 GeV. Remember, that, as well known, 
the wave functions of states bound in the coulomb po- 
tential exponentially decline versus the distance with the 
damping lengths of a n, where n is the principal quan- 
tum number, so that for the bag-charmonium this length 
is less than 1 fm at n — 1, 2. Supposing that the bag it- 
self has the size somewhat greater than 1 fm, we draw the 
conclusion that the approximation introduced for the po- 
tential as the sum of coulomb term and bag mass is quite 
applicable for the low-lying levels of bag-charmonium. 



B. 



States at n = 1, 2 



Fixing the model parameters by means of considering 
the splitting between the two lowest vector and pseu- 
dovector levels of bag-charmonium at n — 1,2 allows us 
to make spectroscopic predictions for further members of 
family with the same values of principal quantum num- 
ber. The uncertainty of such the prediction is about 5 
MeV. The state masses are listen in Tab. [TTJ Emphasize 
that for the coulomb- like states the splitting of levels due 
to the forces depending on the quark spins approximately 
decreases with the growth of principal quantum number 
by the scaling law 1/n 3 , that reasonably agrees with the 
data in Tab. M 





n = 1 


n = 2 


o-+ 


3678 


3865 


\— 


3770 


3876 


2++ 




3875 


1++ 




3872 


0++ 




3868 


1+- 




3873 



TABLE II: Masses of bag-charmonium at n — 1, 2 in MeV. 

Two of listened states have been associated with the 
following ones: 

1~ [l 3 Si] i — * V(3770), l ++ [2 3 Pi] i ► X(3872), 

the quantum numbers J of which are established ex- 
perimentally while the exotic level X(3875) with uniden- 
tified quantum numbers could get the two following pre- 
scriptions: 



l—^Sx], 
2++[2 3 P 2 ]. 
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C. States at n > 2 

The accepted approximation for the potential as the 
sum of coulomb attraction and bag mass is valid at 
n = 1,2, but it can be modified in the case of higher ex- 
citations because of perturbation responsible for a "soft" 
transition into the regime of confinement, since at n = 3 
the system composed of charmed quark and antiquark 
has got the size about 1 fm, i.e. it is close to the confine- 
ment scale, and hence, to the bag size itself. Therefore, 
we have to account for some other nonperturbative terms. 

So, the presence of gluon condensate, i.e. an external 
field, leads to the contribution caused by the second order 
of perturbation theory due to the chromoelectric dipole 
interaction to the next-to-leading order in the velocity of 
heavy quarks, that has the form 7 



5M n = n n (an) 2 



(27) 



at 



.4 



(a s G 2 ) 
5E n ■ 



where (a s G 2 ) ~ Aq CD is the gluon condensate, SE n ~ 
TOQ(|a!g ff ) 2 /4n 2 is the characteristic energy of binding 
the heavy quarks, and A denotes the numerical dimen- 
sionless factor caused by color and spatial effects. Then, 
at Aqcd ~ 0.2 GeV in the case of bag-charmonium 



7 We deal with the general formula for the correction to the en- 
ergy SM ~ (V) 2 /SE with the perturbation due to the electric 
dipole V ~ g s r£, so that the charge squared gives the coupling 
constant g 2 i— > a s , the square of chromoelectric field is reduced 
to the square of gluon stress E 2 i— > — G 2 ^, and the characteris- 
tic size gives the Bohr radius of n-th excitation r t— > an. Ac- 
tually, the strict calculation of nonperturbative term due to the 
gluon condensate as well as the complete analysis of nonrelativis- 
tic quark-antiquark system with coulomb interaction was done 
in |28l . |29| . l30t |31H for the charmonium and bottmonium in the 
framework of multipole expansion of QCD [lj. The authors used 
exact formulae for the Green function of color-octet state, so that 
8M n i = n 6 mQ{a s G 2 ) a 4 7re n i/16. At first, this correction gives 
sixth power of principal quantum number, in fact, because of 
averaging (r) >— » an 2 valid for the coulomb potential (the exact 
formula (r) = a[3n 2 — 1(1 + l)]/2 includes the dependence on the 
orbital momentum). Second, the numerical dimensional factor 
in front of correction is certainly definite in such the estimate. 
Third, one gets the opportunity to take into account for the or- 
bital quantum number in terms of factor e n ; ~ 1.5, which value 
is also strictly known. However, wc do not transfer such the 
analysis to the case of considering the bag-quarkonium, because 
we suggest that the bag essentially distorts the propagation of 
color-octet state, since it cuts off the wave functions in infrared, 
especially in the case of subtraction between the quark and anti- 
quark. In this respect, our approach to the problem is less strict, 
but, to our opinion, it is more close to the actual physical situ- 
ation, so that the fourth power in the dependence of correction 
on the principal quantum number is more realistic, while the nu- 
merical dimensional factor is fitted phenomenologically, because 
the modification of wave functions or Green function in infrared 
can essentially "renormalize" this factor. 



we find the estimate K n ~ An 2 mQa 2 (a s G 2 ) ~ 
An 2 xlO -3 GeV 3 . Remember that the positive value 
of gluon condensate corresponds to "negative" square 
of chromoelectric field; the intermediate state of quark- 
antiquark pair is the color octet with the subtraction be- 
tween the quarks, hence, it has the energy greater than 
the bound states under study, that leads to positive value 
of correction to the mass of state 8 . In addition, the ef- 
fective constant related with the gluon condensate can 
depend on the fact that the chromoelectric field is par- 
tially screened by the bag, which is the infrared object 
itself. This fact can lead to a dependence of effective 
value related with the condensate, that has the step-like 
form changing the value at a critical distance given by 
the bag size, so that k i— » n(l + £'0[n i , — n]). 

Let us evaluate the effect of perturbation introduced 
in the case of neglecting the spin-dependent forces 9 by 
setting 

k = kq(1 + — n]J, 

n n =n 2 K, KQ = l.lxlCT 4 GeV 3 , (29) 

e = 0.3, 4 < < 5. 

The value of kq is in fact fixed by the position of higher 
excitations. It is important that the functional depen- 
dence of correction on the principal quantum number 
truly reflects the structure of extraordinary states in the 
charmonium family. The introduction of nonzero value 
for the parameter e allows us to slightly displace the po- 
sitions of levels at n = 3, 4 by 10 and 30 MeV, respec- 
tively, that is comparable with the accuracy calculation 
of methods generally involved in the framework of poten- 
tial models. Nevertheless, the physical meaning reflected 
by parameters £ and n* rs 4.5 is quite clear: the degree 
of screening the chromoelectric field by the bag is of the 
order of 30 %, and the bag size is about nti ( « 1.5 — 1.6 
fm. 



8 Correction II27I I can be roughly described by introducing the per- 
turbation potential V = n'r 2 , so that for the coulomb functions 
of initial states the shift of energy is equal to 



(V) 



- a 2 n 2 [5n 2 + 1 - 3L(L + 1) 



(28) 



which repeats the general behavior versus the principal quantum 
number, of course. Eq. H28H supposes that the correction to the 
energy is determined by one and the same value of parameter k' 
independent of the state, that is certainly too strong suggestion. 
Though, the final expression can be treated as the factorization 
of overall factor, while the residual dependence versus the pa- 
rameters of state is close to the fourth power of principal quan- 
tum number. In H28H one could account for correction variation 
because of nonzero value of orbital momentum, but this piece 
of subtlety is evidently beyond the accuracy suggested in the 
derivation of correction itself. 
9 As we have already mentioned above, the splitting of higher exci- 
tations is significantly suppressed by the factor of 1/n 3 , so that 
the variation of estimates due to the spin-dependent forces at 
n = 3 is less than 3 MeV, that is certainly below the uncertainty 
in the formula for the interpolation used. 
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Then, we find that at n = 1,2 the contribution of per- 
turbation is less than 1 MeV, while for the higher exci- 
tations we get predictions corresponding to exotic states 
of charmonium: 



M bag n " 31 = 3939 MeV 

M b^ n=41 = 4037 MeV 
M^ n=51 = 4156 MeV 
M^ =6] = 4424 MeV 



V(3940), V(3940), 
Z(3940), 

^(4040), 

^(4160), X(4160), 
^(4415). 



(30) 



Hypothetic higher excitations can suffer from the influ- 
ence of linearly raising term in the potential, so that their 
masses could exceed the naive estimate at n > 5, say, for 

M blg =6] ~ 4 - 42 GeV ( thc 

size of quar-antiquark sys- 
tem exceeds 2 fm) , but such the high values of excitation 
energy signal that thresholds of some other effects could 
open. Particularly, the excitation of bag could take place, 
that we discuss in the next subsection. 



D. Exciting the bag 

The method of estimating the mass of basic state for 
the bag in the framework of constituent model for ir- 
and p-mesons suggests that the bag could have excita- 
tions. The first excitation is related with the system of 
a-mesons. The spin-average state has the mass 10 equal 
to 



m(lP) = i (5m a2 + 3m ai + m ao ) » 1260 MeV. (31) 



Therefore, the mass of low-lying excitation of bag is equal 



to 1 



E' hag * 1105 MeV. 



(33) 



Then, we easily determine the increment of repeating the 
low-lying states in the system of bag-charmonium 



6M hag = E' hag - P bag a 490 MeV, 



(34) 



so that we get the prediction permitting the direct ver- 
ification by the comparison with known exotic states of 
charmonium family, 



MgJl 3 Si]' a 4260 MeV i-> V(4260), 
MgJ&Si]' a 4366 MeV i-> V(4360). 



(35) 



The next excitation of bag is the straightforward ana- 
log by the analysis of tt and p system, namely, by con- 
sidering the radial excitation 25: 7r(1300) and p(1450). 
Then, in accordance with the experimental data 12 



Ka g ~ 142 ± 25 MeV. 



(36) 



Therefore, the masses of vector states in such the bag are 
equal to 



M£l g [l 3 Si]" « 4575 ± 25 MeV 



M«L[2 3 5i]" a 4680 ± 25 MeV i-> y(4660). 



(37) 



Thus, the model of bag-charmonium allows us to to- 
tally describe the wide spectrum of experimental data 
on the spectroscopy of extraordinary states in the char- 
monium family by following the qualitatively clear as- 
sumptions and fixing the parameters from the measured 
splitting between the vector and pseudovector states. 



However, in contrast to the 5-wave state, the presence 
of P-wave suggests that the constituent quarks get the 
kinetic energy of orbital rotation, 



E n 



bit 



2m le dr 2 ' 

where the reduced mass is equal to 

^rcd 2 ^const ~ 4 Pbag- 

At L = 1 and L ~ m rc( j r, we find 

Porbit ~ m lc d « - P b ag ~ 155 MeV. 



(32) 



10 We use the prescription of n r l, where n r is the radial quantum 
number. 



E. The spectrum of charmonium family 

Thc final spectrum of states for the bag-charmonium is 
shown in Fig. [5] in comparison with the experimentally 
measured positions of levels in the charmonium family 
including the extra states arranged by the total momen- 
tum, spatial and charged parities, J PC , if the quantum 
numbers are empirically established or in the case of bag- 
charmonium with a sizable splittings depending on the 
spin, i.e. at n = 1,2. The higher excitations at n > 2 



One could get the similar estimate, if one starts from the calcu- 
lation of mass for the P-wave state of D-meson, so that in the 
same approximation we obtain mrj (IP) = m c + i _B^ a + E^ b ., , 

where P D rbit « rnP, d = m c m CO nst/(m c + m CO nst) « 260 MeV, 
while the experimental data are not complete in order to cal- 
culate thc position of spin-average level, but we could hold the 
trend by setting m™ p (lP) (5m D . + 3m Dl )/8 « 2445 MeV. 
12 The mass of 7r(1300) is measured with uncertainty of 100 MeV. 
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E, MeV . 



4600 





r<4<m_j. _ 

i Hi 






ill [441 5) 6S 
Y{4366) 2S> 

Y(4260) IS , 






(If (4160) ss 




D-JY 


<fi(4040) 4S 




pja, 


3S 


mom Jre*^ 39 *" 




_25 


A'(_*S72) /P X^JfiZJ) 


DD 








<lr(2S) 








XcUFi 















0>- 1— f- / 4+ 2 ++ ? J r + J re 



FIG. 5: The spectrum of charmonium family arranged by the quantum numbers J : solid lines mark the experimental data 
[3^ ]; dashed lines show the predictions in the model of bag-charmonium (see details in the text); dotted lines are positioned 
at the thresholds for the production of mesons with the open charm. The shaded bands give uncertainties in the case of both 
experimental data and theoretical predictions. 



are marked by the dashed line drawn in the limits from 
1 to 2 ++ in order to clearly point to the fact that such 
the prediction concerns for the permitted higher values 
of orbital momentum, not only the vector states. In that 
case one has to take into attention that the 5-wave label 
nS refers to the pseudoscalar and vector mesons, only. 

The levels of bag-charmonium are marked by the ra- 
dial and orbital quantum numbers: n r l. The prime in the 
notation means that the bag has got its first P-wave ex- 
citation, and the double prime points to the first S-wave 
excitation of bag. 

In favor of state identification accepted, there are the 
following facts: 

• the anomaly in the distribution over the invariant 
mass of two pions in the decay -0(3770) — > J/tpmr, 



• the preference (or the evident selectivity) in 
two-pion transitions of F(4260) — > J/xjjim and 
F(4360) -> iP'tttt for the IS- and 25-states of char- 
monium, respectively, 

• the analogous selectivity in the two-pion transition 
of F(4660) — > ij/irir for the certain 2S'-level of char- 
monium. 



III. TWO-PION TRANSITIONS 

The interaction of compact color-singlets composed of 
heavy quark and antiquark, with "soft degrees of free- 
dom" , for instance, with the low energy pions, is con- 
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structed in the framework of multipole expansion of QCD 
[l[ , so that in the leading order over the heavy quark ve- 
locity the transition of bag-quarkonium to the state of 
quarkonium with the emission of two pions posed in the 
isospin-singlet and S'-wave state, takes place due to the 
second order over the chromoelectric dipole moment with 
the amplitude 13 ^4[bag-(QQ) 1 — > (QQ) 2 7T7r] = A equal to 

A = ( 7 r 7 r(QQ) 2 |« ?s 2 q4r 7 T a ar b r /3 |bag-(QQ) 1 ), (38) 

where T° is the "difference of generators for the color 
charges" of quark and antiquark: T a = |(*q — m 
terms of Gell-Mann matrices A a = 2t a , Q denotes the 
Green function of color-octet state of quark and anti- 
quark, while subscripts of quark-antiquark state mark the 
set of its quantum numbers. In matrix element (|38p one 
could isolate two following contributions: the first cor- 
responds to the factorization of soft degrees of freedom 



^[bag-CQQ)! -» (QQ) 2 7T7r] 



SF 



Asf, that takes place in 



the leading order over the velocity of heavy quarks, 



A SF = (7T7r|^f»^|bag)x 

((QQh\r 7 T a gT b r \(QQ) 1 ), 



(39) 



so that one could introduce the chromoelectric polariz- 
ability in the transition between the color-singlet states 
of heavy quark and antiquark 



<4 2) - 7 ((QQ) 2 1 r^T a gT a rp \(QQ)i), 



1 



which gives 



A 



SF 



2^ 



(12) /w| ffs 2 £ 7 ^|bag) 



(40) 



(41) 



The second contribution appears due to account for the 
bag annihilation ^4[bag-(QQ) 1 — > (QQ) 2 7T7r] = Aba, 

BA 

which becomes possible in higher order over the heavy 
quark velovity 14 , 



Aba = (ttttI^^IO): 



((QQ) 2 |r 7 T a gT b r /3 |bag-(QQ) 1 ) 



(42) 



that is also reduced to the introduction of analog to the 
chromoelectric polarizability 



7/3 



1 



a 



((QQ) a | r 7 T a QT a r fj |bag-(QQ) 1 ), (43) 



which can be nonzero beginning from the second order 
over the heavy quark velocity, since in the first order 



over the velocity the operators of interaction (dipole mo- 
ments) take nonzero color charge, a ~ 0(vq). The re- 
striction to the .S'-wave vector states of quark and anti- 
quark with polarization vectors ei j2 immediately gives 



a lj3 = aSj/s (ei • e 2 ), 



hence, 



A = ia( 12 )(w| ff 2 ^|bag) 



+ ia( 12 )(w| ff 2 ^|0). 



(44) 



The standard technique 33], presented, for instance, in 
review [17| . allows us to relate the operator quadratic in 
the chromoelectric field with the anomaly in the trace of 
energ y-m omentum tensor in QCD, so that in the chiral 
limit [3J| we get 



(7r7 r| 5 2 ^|0) 



4^ 

' AT 



(45) 



where q = pi + p 2 is the sum of 4-momenta of pions 
in the final state, and q 2 denotes their invariant mass 
squared. Analogous speculations in combination with an 
expansion of "soft" nonstationary state of bag over the 
pion states with positive spatial and charge parities and 
zero charges with respect to gauge interactions lead to 
the conclusion that in the limit of soft chiral pions the 
dominant contribution in the matrix element is given by 
the two-pion projection of bag state and, hence, four-pion 
vertex, i.e. the constant integrated over the parameters 
of expansion for the bag, that is approximated by 



1 



(Trvrl^^lbag 



4^ 

' AT 



(46) 



where 8 defines the complex phase of matrix element (|46[) 
with respect to f|45[) . while the dimensional parameter 
is determined by the bag mass fio ~ -E-bag, i-e. it has 
the magnitude of the order of several Aqcd- Thus, the 
two-pion transitions between the vector S'-wave states of 
quarkonium in the presence of bag in the initial state are 
described by the functional dependence 



A[bag-(QQ) X -> (QQ) 2 7T7r] 



? 2 +M 2 e i5 , 



(47) 



wherein the bag causes the modification by introducing 
the term with /i =/= 0. 

It is well known that the parametrization 15 of (|4"7| 
is quite accurately describes the distribution of pions 
over their invariant mass in the analogous transition 
T(3S f ) — > Y7T7T. At first, this fact signals that T(3S') can 



See the methodic presentation in review fl7l| ■ 
In addition, the propagation of color-octet state of quarkonium 
with subtractive forces suggests that the wave functions of quarks 
are essentially overlapped with the bag, that can result in the 
annihilation of bag. 



The standard parametrization includes the modification allowing 
us to take into account the fine effect caused by the rescattcring 
of pions in the final state [35ll . that produces corrections, which 
accuracy is graded by the empirical uncertainty of real data on 
the spectrum of distribution over the invariant mass of pions ^Qj . 
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be treated as the basic vector state of bag-bottomonium. 
Second, it shows that fitting the parameters in the second 
term of (|47j) certainly allows us to describe the anomalous 
distribution in the transition of tp(3770) — > J/ipTnr, too. 
In this way, the value of /x is close to 0.7 GeV in the case 
of T(35'), that is consistent with the prediction in the 
model of bag-quarkonium. On the other hand, the ap- 
pearance of such the scale in the transitions between the 
ordinary states of heavy quarkonium with the emission 
of two pions seems to be rather problematic. Though, 
at present, there are various versions for the explanation 
of anomaly in question, which are not reduced to our 
suggestion, of course (see, for example, ^7\). 

Thus, the anomalous distribution over the invariant 
mass of two pions in the transition of bag-quarkonium 
into the heavy quarkonium is consistent with the ac- 
cepted identification of basic vector states of bag- 
charmonium and bag-bottomonium. 



IV. LEPTONIC CONSTANTS 

The leptonic constants of heavy quarkonium can be 
calculated in the potential model by applying the effec- 
tive lagrangian of nonrelativistic QCD (NRQCD), that 
was studied in 38] in the case of heavy quark and anti- 
quark of the same flavor and in [39] ] for the quarkonium 
composed of heavy quark and antiquark of different fla- 
vors. So, the vector currents of quarks Q in nonrelativis- 
tic QCD are related with the currents in the full theory 
by the formula 



J? CD = Qi,Q, 



j- NRQCD = x](jU 



where Q is the field of relativistic quark, while x an d 
<f> demote the nonrelativistic Pauli spinors for the quark 
and antiquark, a^; = a l/ — v l/ (a-v), and v is the 4- velocity 
of quarkonium, so that 



jQCD _ /C(^hard) Mfact) ' J^ R( ^ CD '(^f act ). 



(48) 



Here /Lthard fixes the point of matching NRQCD to QCD, 
Mfact denotes the scale of perturbative calculations in 
NRQCD. 

In the case of quarks of the same flavor, the Wilson 
coefficient /C is known to the 2-loop accuracy SH, S3, 
El 



MMhardi A t factJ= 1 — -- 



3 7T 
a^^hard) 



(49) 



C2(A*hard; Mfact)) 



C2 is explicitly given in [4l|, |42| . The anomalous dimen- 
sion Of /CdUhaxd; /"fact) 



in two loops is given by the expression 16 
7[i] = 0, 

7[2] 



= -Wtt^Cf -C f 



(51) 
(52) 



Initial data for the evolution of /C(/ihard; Mfact) versus the 
scale are determined by the matching conditions fixed at 

M = Mhard EE S3. 

The leptonic constant of vector state with the polar- 
ization A and polarization vector in full QCD 

<0|J, QCZ5 |QQ,AH^/ QQ Af QQ , (53) 

is related with the matrix element of current jNRQCD 

^jNRQCD^^ A) = A{fx) £ X f NRQCD M ^ (M) 

where the potential model of quarkonium with the wave 
function 'I'QQ(r) gives 



NRQCD 



I*qq(°)I- 



(55) 



while the renormalization group factor A(n) equals unit 
at the scale of /i = /j,q, whereat the wave function is 
determined. Since 

IqQ = Jqq D -4(Mfact) • £(Mhard;^fact), (56) 

the anomalous dimension of *4(^f act ) compensates the 
anomalous dimension of /C(/Xhard; Mfact), hence, 



din /it 



-7[2] 



4tt 



(57) 



therefore, 



A{n) = A(no) 



0o + A 



4tt 



47T 



2/3i 



(58) 



The calculation of leptonic constants in wide limits for 
the scales of matching and factorization allows us to find 
the region, wherein the result is stable, i.e. it slightly 
depends on small variations of scales. Then, is set to 
the point of stability. Numerical estimates are consistent 
with the experimental data for the basic vector states of 
charmonium and bottomonium [381 ] . while the leptonic 
constants of excited states in such the approach are sim- 
ply determined from the ratio of wave functions at the 
origin for the given state to the wave function of basic 



= 1^ 7[fc] ( ) > ( 5 °) 



fc=i 



47T 



16 In terms of ordinary notations for the representations of SU(N C ) 
group: C F - -r^- 



N Z — 1 1 

v 1 -,C A =N C> T F = L 
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state due to appropriate rescaling of leptonic constant for 
the basic state in accordance with (|56p . since the renor- 
malization group factors are universal , i.e. they do not 
depend on the excitation number. Then, the problem 
is reduced to the safe calculation of wave functions of 
quarkonium at the origin. Such the calculation involves 
rather a large uncertainty because of modelling the po- 
tential and prescribing the pole masses of quarks depend- 
ing on the model. Thus, one arrives to a sizable methodic 
uncertainty of theoretical predictions. Nevertheless, one 
can reliably state that the low value of leptonic constant 
of "0(3770) is caused by the fact that this state cannot be 
assigned to the S-wave state of charmonium. Therefore, 
it is the exotic state whether it is the result of mixing the 
S'-wave level with the D-wave one, as usually accepted, 
in spite of the fact that the splitting between such the 
nearest levels is about 100 MeV, and hence, the matrix 
element of mixing should get the same order of magni- 
tude, that is quite a large value. 

However, such the framework is not applicable to the 
states of bag-quarkonium, since in this case the joint an- 
nihilation of quark-antiquark pair and bag should take 
place. Then, the knowledge of wave function for the 
quark-antiquark system is not sufficient. In addition, the 
wave function could be essentially modified during the 
annihilation with the bag. So, let us explore the method 
of quasilocal sum rules of QCD [2(J. The basic idea of 
such the approach is reduced to the following: 

• Consider the transversal part of correlator for the 
vector currents of heavy quarks, namely, the real 
part of correlator at the invariant mass set to zero 
as well as the derivatives of k-th order, i.e the mo- 
ments of real part at zero. 

• In the case of heavy quarks, the bound states are 
positioned in the narrow gap which width is sup- 
pressed in comparison with the heavy quark mass 
niQ, so that for the low values of moment numbers 
one can transform the summation over the reso- 
nances to the integration over the density of states 
dn/dMn, where n denotes the number of vector 
state in the direction of mass increase 17 , with ac- 
curacy up to 0(1/ mo). 

• In the leading order over the inverse mass 1 / toq at 
low values of moment number, calculating the cor- 
relator can be restricted by the contribution of non- 
relativistic quarks with account for their coulomb- 
like interaction. 

Then, at the points of resonances the local equality of 
theoretical correlator to the empirical one approximated 
due to the introduction of density of states leads to the 



relation for the leptonic constant of n-th state fa 

Jl^ c dMn 
Mfi dn ' 

where the constant is defined by the formula 



/ 2m c 



7T \ Ma 



^sys i 



(59) 



(60) 



with an effective constant of coulomb exchange 18 a c s , the 
factor of loop corrections K. taken to the 1-loop accuracy 
from (|49p and a systematic correction generally depend- 
ing on the moment number 



Z 



NR 
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introduced in the both nonrelativistic approximation of 
quark correlator and density of states for the integral 
representation for the sum over the resonances. In the 
limit of heavy quarks, i.e. by neglecting logarithmic and 
power corrections of the form ln/x/mQ and I/ttiq, the 
effective constant C is independent of quark flavor. It is 
spectacular, that empirically for the finite masses of real 
c- and 6-quarks and basic vector states, we have got 19 

so that the scaling relation of (fS"9"l) takes place. 

If the heavy quark and antiquark states are certainly 
the levels of nonrelativistic heavy quarkonium, then the 
density of states is, in practice, universal, since the po- 
tential is close to the logarithmic dependence, so that 
independently of the heavy quark flavor we get 



dM h 2T 



dri 



(61) 



17 We insert the tilde symbol in order to distinguish the number 
from the principal quantum number n used above. 



The effective constant of coulomb exchange as well as the quark 
mass in this approximation do not necessary coincide with their 
analogs in the framework of potential approach. Though, the 
mentioned differences become less with the increase of quark 
mass, so that for the bottomonium these constants coincide with 
the accuracy below 10%, indeed, if one does not take into ac- 
count for the difference in the values of 6-quark pole mass, while 
the account for the renormalization group evolution of coupling 
constant versus the characteristic momentum inside the quarko- 
nium, hence, versus the quark mass, leads to that the effective 
constants actually coincide with each other with the accuracy 
better than 2%, that is essentially underestimate the accuracy of 
assumptions made in the way of deriving these estimates. In ad- 
dition, the strict analysis of charmonium and bottomonium made 
in the framework of assumption on the dominant coulomb inter- 
action with account for the perturbative and nonperturbative 
corrections in [28, 29, 30, 31] led to the values of pole masses of 
charmed and bottom quarks, which are consistent with the values 
obtained in the present paper within the accuracy of estimates. 
The effective constants of coulomb interaction take similar val- 
ues, too. 

The exact value can depend on the choice of quark mass, that 
leads also to a variation of systematic factor caused by the non- 
relativistic approximation in the calculation of coulomb contri- 
bution near the threshold. 
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where T is the average kinetic energy, i.e. the main 
parameter of logarithmic potential. Therefore, for the 
heavy quarkonium the scaling law takes place 



const. 



(62) 



valid for J/0 and T. By the way, since the width of 
vector quarkonium decay r[lqq — > e + e~] is related with 
the electric charge of quark eQ and its leptonic constant 
by the expression 



47T 

~9~ 



f 2 
J n 

M ft ' 



for the standard states of heavy quarkonium with the 
regular spectrum under condition (|6ip one could expect 
the universal law 



-5- Tnfi — const. 



(63) 



Relation ([55)) is also applicable, if there are exotic, ex- 
traordinary states with the hidden charm, too. Its phys- 
ical meaning is simple: the contributions of resonances 
into the correlator are actually fixed by the dislocation 
in accordance with the observed density of states. There- 
fore, it is enough to carry out the analysis of data on the 
spectrum of vector states with the hidden charm and, 
then, to test the direct consequence of (|59|) in the form 
of 



1 „ dn 

— Ifi — — = const. 

en dMn 



(64) 



The accuracy of analyzing the spectrum of vector 
states is essentially restricted by the fact that the lev- 
els are discrete. The result of cubic interpolation for the 
mass spectrum with the further calculation of dMn/dn 
versus n being the number of vector state as given in Fig. 
[5] by the mass increase, is shown in Fig. [5] 

Note that the interpolation of spectrum itself versus 
the resonance number leads to discontinuities of the mass 
derivative with respect to the number exactly at the 
points, which refer to the values of state density in in- 
terest. Therefore, we have made the cubic interpolation 
of the state density itself once more with different in- 
crements in the number in order to smoothen the dis- 
continuities 20 and show uncertainties of estimating the 



20 The artefact of function discontinuity points to the fact that the 
derivative of spectrum with respect to the number of resonance 
permits a displacement of point on the number axis in the limits 
Sn = 0.5, so that we have also applied the method of smoothing 
by averaging the values of spectrum density taken left and right 
in vicinity of resonance, that agrees with the uncertainties of es- 
timation itself. Then, the smoothing is slowly modify the values 
of spectrum density in the points, where it fluctuates weakly, i.e. 
the spectrum is rather regular, and in addition, it allows us to 
draw the limits of approach uncertainties. 
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FIG. 6: The result of cubic interpolation of vector state spec- 
trum of charmonium family in accordance with predictions 
of model with the consequent calculation of AMn /dn (dashed 
line) in comparison with the smoothed cubic interpolation of 
the state density involving the increment steps equal to 1.5 
and 0.8 in h (solid and dotted lines, correspondingly). 



quantity of dM^/dn. In figure we see that the most re- 
liable results are given at n = 6,7 (states -0(4040) and 
■0(4156)), whereat the spectrum is the most regular in 
the sense of its smoothness by incrementing between the 
nearest states. The uncertainty grows at h = 3 (the state 
0(3770)), whereat the mass increment is extremely inho- 
mogeneous, as well as at the end of interpolation interval, 
i.e. at the important positions of states with n = 1,2, 
mesons J/ip and 0(25). The absolutely analogous situ- 
ation takes place at other choice of interpolation power, 
namely, at linear and quadratic interpolations, for in- 
stance. This fact means that the investigation of relation 
(f6"4"| should be made by fixing the constant at h = 6, say, 
that we will perform further. In this way, an essential un- 
certainty in the calculation of leptonic constant of basic 
state J/ijj appears. In order to decrease the uncertainty 
of interpolating towards J/0, we will take into account 
that this state is certainly the heavy quarkonim, hence, 
the state density in vicinity of n = 1 is rather close to the 
value predicted in the potential model. This fact can be 
naturally involved by introducing an auxiliary state with 
h = and mass M « Mj /lP - 2T « 2296 MeV according 
to the linear extrapolation, that we actually make. 

Thus, relation (f6~4"]) leads to estimation of leptonic con- 
stants shown in Fig. [71 We see that the scaling rela- 
tion quite successfully describes the behavior of leptonic 
widths for the observed states. Since the relation is uni- 
versal, and it is valid irrespectively of the role played by 
infrared phenomena such as the bag, there is no need to 
prescribe the state of 0(3770) as the ID level of char- 
monium with the appropriate mixing with the 2S'-level 
in order to produce nonzero value of leptonic constant 
(the wave function at the origin). Moreover, the small 
difference between the leptonic constant of 0(4040) and 
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F#, keV 




J/iii $(25) ${3770) 2S 3S ${4040) ${4156) Itf 2S' ${4415) 



FIG. 7: The leptonic constants of vector states in the family of charmonium T„ in keV versus the resonance number h (solid 
line). The shaded region marks the systematic uncertainty of predictions, which is caused by the interpolation of state density. 
The dots denote the experimental data with the uncertainties shown as vertical lines. 



'0(4156) looks quite exotic, if one treats them as the S- 
wave levels of charmonium with the respective numbers 
h i— ► 3, 4, while the presence of extra states due to the 
bag-charmonium leads to the both equalizing the ratio 
of leptonic widths and more natural values of absolute 
values: T§ w Similar statements can be written as 

concerns for the state of ^(4415). 

We do not present a complete description for the pro- 
cedure of interpolating the spectrum density, since Fig. 
[7] clearly shows that the method inherently involves irre- 
ducible uncertainties, which do not permit to make the 
theoretical predictions with the higher precision 21 , com- 
parable with the accuracy of experiment. Nevertheless, 
the exploration of quasilocal sum rules by introducing the 
density of bound states with the hidden charm allows us 
to establish the tendency in the character of dependence 
of leptonic constants on the number of excitation over 
the basic state. 

Thus, the extraordinary states in the family of char- 
monium as dictated by the model of bag-charmonium, 
do not contradict with the analysis of know leptonic con- 
stants of vector states, and moreover, they permit to im- 



21 For instance, we could quite smoothly fit the spectrum of vector 
states with the hidden charm by means of special set of basic 
functions, so that the spectrum density will repeat the exper- 
imental dependence of leptonic constants versus the resonance 
number. However, the procedure of selecting the basic functions 
is arbitrary itself, so it involves almost uncontrolled methodic un- 
certainty, restricted by the constraint of acceptable smoothness 
of the fit, only. 



prove the systematization of those constants (or leptonic 
widths). 



V. CONCLUSION 

In the present paper we have formulated the poten- 
tial model of bag-charmonium, which has allowed us to 
successfully identify extraordinary states in the family of 
charmonium in terms of spectroscopy, to substantiate the 
anomaly in the distribution over the invariant mass of two 
pions in hadronic transitions such as ^(3770) — > J/i/)ittt 
and T(3S') — * T7T7T by assigning the initial states as the 
vector bag-quarkonia and to establish the systematic reg- 
ularity of leptonic constants for the vector states in the 
charmonium family due to the levels of bag-quarkonium. 
The regularity has got the form of scaling law derived 
from the quasilocal sum rules of QCD. 

However, such the model predicts some new states yet 
not observed empirically. 

As concerns for the mechanisms of production of ex- 
otic states, one could stress the following: in decays of 
heavy hadrons with the light valence quark, the produc- 
tion of bag-charmonium goes almost the same way as the 
production of charmonium, namely, at rather short dis- 
tances about the charmonium size, the compact state of 
charmed quark and antiquark is formed, while the bag 
itself already exists due to the presence of light valence 
quark in the initial state, so that the bag introduces 
an appropriate form factor, describing a possible tran- 
sition of bag from the rest in the initial state to a motion 
joint to the quark-antiquark pair with the hidden charm. 
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An analogous form factor, probably, should appear at 
the production of bag-charmonium in hadronic collisions, 
whereas the order of magnitude for the cross section is 
close to the production cross section of charmonium. In 
leptonic collisions, the main role is played by the leptonic 
constants studied in Section ITVl 

In the present paper we, in fact, have not discussed 
any channels of decays for the bag-quarkonium except 
pointing to the selectivity in the two-pion transitions of 
bag-charmonium to the charmonium states with identical 
quantum numbers of quark- antiquark pair. 

In this respect, it is appropriate once more to refer 
to review [19(, wherein there is a comprehensive bibli- 
ograthy of original articles concerning for the models of 
exotic quarkonium states, their decays and production 
modes as well as a comparative description of various 
mechanisms of forming the extraordinary states in the 
charmonium family. We point to some pioneer papers, 
wherein authors introduced such the notions into the sci- 
entific language in the field of quark dynamics of exotic 



states with the hidden charm as follows: 

• the hadronic molecule 0, |45[ or deuson (in analog 
to the deutron) [46| , composed of charmed mesons 
with the positive and negative charm DD, 

• the hybrid pTl |48| containing the quark-antiquark 
pair in the color-octet, which is coupled with a com- 
pact gluonic lump to form the colorless state, 

• the tetraquark [49j composed of diquark and dian- 
tiquark. 

In addition, it worth to mention the recent paper on 
the systematization of heavy quarkonium states - char- 
monium and bottomonium - in the scheme of Regge 
trajectories [50| . Studying various color structures of 
tetraquark has been given in [5lj . 

This paper is partially supported by RosAtom, con- 
tract # N.4d.47.03.08.069, and RFBR, grant # 07-02- 
00417. 
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